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Abstract The collisional rates associated with the isotropic 
velocity moments (V^' ) and the anisotropic moments {V^'Yi) 
and {V^''{ViVj — d^^V^dij)) are exactly derived in the case 
of the inelastic Maxwell model as functions of the exponent 
r, the coefficient of restitution a, and the dimensionality d. 
The results are applied to the evolution of the moments in 
the homogeneous free cooling state. It is found that, at a 
given value of a, not only the isotropic moments of a degree 
higher than a certain value diverge but also the anisotropic 
moments do. This implies that, while the scaled distribu- 
tion function has been proven in the literature to converge to 
the isotropic self-similar solution in well-defined mathemat- 
ical terms, nonzero initial anisotropic moments do not decay 
with time. On the other hand, our results show that the ratio 
between an anisotropic moment and the isotropic moment 
of the same degree tends to zero. 

Keywords Inelastic Maxwell model • Collisional rates • 
Homogeneous cooUng state 



1 Introduction 

The prototypical model of a granular gas consists of a sys- 
tem of (smooth) inelastic hard spheres (IHS) with a con- 
stant coefficient of normal restitution < a < 1 fTSi . Under 
low-density conditions, the one-particle velocity distribution 
function /(r, v;f ) obeys the (inelastic) Boltzmann equation. 
On the other hand, because of the intricacy of the colli- 
sion operator, one has to resort to approximate or numeri- 
cal methods to get explicit results, even in the elastic case 
(a = 1). The main mathematical difficulty lies in the fact 
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that the collision frequency of IHS is proportional to the rel- 
ative velocity of the two colliding particles. As in the elas- 
tic case II20II29I . a significant way of overcoming the above 
problem is to apply a mean-field approach whereby the colli- 
sion frequency is replaced by an effective quantity indepen- 
dent of the relative velocity. This defines the so-called in- 
elastic Maxwell model (IMM), which has received much at- 
tention in the last few years, especially in the applied math- 
ematics literature (see, for instance, l.2..3..4^6^7..8. .10lfTTlfT2l 
fl31fT4lfl7in^l2T]l221l23ll24ll26ll27ll and the review papers ^ 

mmmEa). 

Although the Boltzmann equation for the IMM keeps be- 
ing a mathematically involved nonlinear integro-differential 
equation, a number of exact results can still be obtained. In 
particular, the collisional velocity moments of a certain de- 
gree k can be exactly expressed as a bilinear combination 
of velocity moments of degrees k' < k and k" = k — k' . Of 
course, the terms with k' = k or: k" ~ k sue products of a 
moment of degree k and a coefficient proportional to density 
(moment of zeroth degree). We will refer to the latter co- 
efficient as a collisional rate. While all the collisional rates 
have been evaluated in the one-dimensional case Q, to the 
best of our knowledge, only the ones related to the isotropic 
moments of any degree Il5l l22l and those related to isotropic 
and anisotropic moments of degree equal to or smaller than 
four 1241 have been obtained for general dimensionality d. 

The aim of this paper is to derive the collisional rates as- 
sociated, not only with the isotropic velocity moments (V"'), 
but also with the anisotropic moments (V^'Vi) and (y^''{yiVj - 
d^^V^dij)). This is done by a method alternative to that fol- 
lowed in Refs. Il5]l22l for the isotropic moments. The knowl- 
edge of the above collisional rates is applied to the study 
of the time evolution of the moments in the homogeneous 
cooling state (HCS). It is known that the isotropic moments, 
scaled with respect to the thermal velocity, diverge in time 
beyond a certain degree that depends on a, as a consequence 
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of the algebraic high- velocity tail exhibited by the HCS self- 
similar solution Il4l l2ni22l . The relevant finding of our study 
is that, at a given value of a, also the anisotropic moments 
diverge beyond a certain degree. This might seem to be a 
paradoxical result in view of the mathematical proofs, both 
in weak ll6ll8ll9lfT0lfTn[T2l and strong |[T7ll23l senses, that 
the scaled distribution function /* tends for long times to- 
ward the isotropic HCS self-similar solution (poo for any ini- 
tial state (isotropic or anisotropic) with finite second-degree 
moments. The solution of the paradox lies in the fact that the 
above convergence properties do not imply that any moment 
of /* of degree higher than two should converge toward 
the corresponding moment of (j)oo. In fact, our results pro- 
vide a counter-example of that strong moment-based con- 
vergence property. On the other hand, we show that the ratio 
between an anisotropic moment and the isotropic moment 
of the same degree goes to zero. 



2 The inelastic Maxwell model 

In the absence of external forces, the inelastic Boltzmann 
equation for a granular gas reads ifTSi 



(a,+vV)/(r,v;r)=7[v|/,/], 



(1) 



where 7[v|/,/] is the Boltzmann collision operator The form 
of the operator J for the IMM can be obtained from the form 
for IHS by replacing the IHS collision frequency (which 
is proportional to the relative velocity of the two colliding 
particles) by an effective velocity-independent collision fre- 
quency Is). With this simplification, the velocity integral 
of the product h{\)J[\\f,f], where h{\) is an arbitrary test 
function ("weak" form of 7), becomes 



dvi/i(vi)7[vi|/,/] = ;^ / dvi / dv2/(vi)/(v2) 



X /da[/i(v'/)-/<vi)], 



where 



Vi 



> 



-a)id-g)a 



(2) 



(3) 



denotes the post-collisional velocity, g = Vi — V2 being the 
relative velocity and a <l being the constant coefficient of 
restitution, n is the number density, 12^ = 27t'''^/r{d/2) is 
the total solid angle in d dimensions, and V is the effective 
collision frequency, which can be seen as a free parameter in 
the model. In particular, in order to get the same expression 
for the cooling rate as the one found for IHS (evaluated in 
the local equiUbrium approximation) the adequate choice is 

HHiia 



d + 2 

V = — r— Vo, Vo = 



4X2,/ 



V^{d + 2) 



no 



d-\ 



(4) 



where CJ is the diameter of the spheres, m is the mass, and T 
is the granular temperature. However, the results derived in 
this paper will be independent of the specific choice of Vq. 

In the case of Maxwell models (both elastic and inelas- 
tic), it is convenient to introduce the Ikenberry polynomi- 
als 1291 F2,|,„,..,- (V) = y2'-y,^,^ ,.^(v) of degree k^2r + s, 
where V = v — u(r) is the peculiar velocity, u(r) being the 
mean flow velocity. The .rth-degree polynomials i'i,,-2...ij(V) 
are obtained by subtracting from V,-, V,^ . . . Vt^ that homoge- 
neous symmetric polynomial of degree s such as to make 

In particular, for s = Q, 1, and 2 one has 



<-iii2...is(y) vanish upon contraction on any pair of indices. 



i2Ho(v) = y2'-, y2.|Kv) = v-'v„ 



>2,-|,7(V) 



V 



2f 



y^Vj-^v^^ij 



(5) 



(6) 



Henceforth we will use the notation M2, i^ and J2r\s^ where 
s = i\i2 ■ ■ ■ is, for the moments and collisional moments, re- 
spectively, associated with the polynomials i2r|.s(^)- Note 
that the collisional moments are defined by Eq. ^ with 

h -^ Y2r\s. 

As said before, the mathematical structure of the Maxwell 
collision operator implies that a collisional moment of de- 
gree k can be expressed in terms of velocity moments of a 
degree less than or equal to k. More specifically. 



J2r\s^ -V2r\sM2r\s+ ^ V,."|.v'.v".v^2,-'|.v'^2r"|.«" 



(7) 



where the dagger in the summation denotes the constraints 

2{r' + r") + i' + / = 2r + s, 2r' + i' > 2, and 2r" + s" > 2. 
Since the first term on the right-hand side of Eq. (|7]) is linear, 
then V2r\s represents the collisional rate associated with the 
polynomial Y2r\s{y)- In particular. 



V2|0 = ^^(1-«')V0, 



(8) 



(l + aMd+l-a) il + af 

V0|2 = ^ \-, ^vo = V210 + ^ . ^ Vo. (9) 

2fl '4 

The quantity V2|o is actually the cooling rate, i.e., the rate of 

change of the granular temperature due to the inelasticity of 

collisions. In general, it is possible to decompose V2, is as 



V2,-|.s 



2r- 



-V2|0 + fi>2,-|. 



(10) 



The first term is the one inherent to the collisional cooling, 
while the second term (02^1.5) can be seen as a shifted colli- 
sional rate associated with the scaled moment 

M2r\s 



Ml 



2r\s 



n(2r/OT)(2'-+.v)/2- 



(11) 



The explicit forms for the collisional rates V2r\s ^nd the 
X coefficients appearing in Eq. (|7]) have been evaluated in 
Ref. II24I for 2r + s <4 and general d. 
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3 Evaluation of V2, O) V2, 1, and V2, 2 

The aim of this section is to evaluate the collisional rates 
V2r|o, V2r|i, and V2,|2 associated with the polynomials (|5]l 
and ^ as functions of the coefficient of restitution and the 
dimensionality. The procedure consists of inserting the poly- 
nomials h = i2,.|o, h = Y2r\i, and h = i2r|i; into Eq. (|2|i and fo- 
cusing only on the term proportional to the moments M2,-\q, 
M2r\i, Sind M2r\ij, respectively. 

Let us describe the method with some detail in the case 
of V2r\o- From the collision rule (O one gets 



(!=1 



vr-'--v^^ = i^rw^''-'\i+aY{d-g) 



r 

l + a 



(a-g)-(a-vi) 



(12) 



//2r 



This equation expresses the difference V[ — Vi as a lin- 
ear combination of terms of order V^' Y^' with ri + r2 = 2r. 
Now, we need to extract those terms of order Vj^'' and V^'' 
only. The terms of order V^'' are obtained from Eq. ( fT2b by 
formally replacing g ^^ Vi, while the terms of order Y^' are 
obtained by formally replacing g ^- — V2 and taking the term 
corresponding to £ = r in the summation. Therefore, 



l + a 



(a-V2)2'- + 42,-|o(Vi,V2),(13) 



where 42r|o (V 1 , V2 ) denotes terms of order W\ ' Vj ' with r\ + 
r2 = 2r, r\ 7^ 0, and r2 7^ 0. When inserting Eq. ( fT3] ) into Eq. 
(|2|i, and ignoring A2r\o^ i^'Sj), we obtain — V2,|oA^2r|o with 
the following expression for V2r|o- 



V2r|0 



V 



t C) '-)'(¥ I ^' 



l + a 



Br 



(14) 



where B( = /dff (ff • g)^'' = In'^'^-^'^l^r (£+5) /F (£ + f) . 
Equation ( fT4l i can be rewritten in a more compact form as 



V2r|0 = Vo- 



(i + 2 



1- 



1 + a 






— 2^1 



1 d 



(15) 



where (0)^ denotes the Pochhammer symbol H], 2^1 {a,b; c;z) 
is the hypergeometric function fl], andz= (l + a)(3 — a)/4. 
Equation ( fTST i agrees with the result derived by Ernst and 
Brito II22I by a different method. 



Proceeding in a similar way, and after lengthy algebra, 
one can evaluate the collisional rates V2,-|i and V2,|2- The re- 
sults are 



V2H1 = vo- 



d + 2 



1- 



1 + a 



2r+I /3\ 
\2>'' 



d{\ + i) 



1 d 



l+a / 3 d+2 
■2F1 ~r,-;^^^;z 



2d 



2' 2 



(16) 



V2,-|2 = Vo- 



d + 2 



l + ax2('-+i) ,+ 1 



il)r 



d{l+d/2){2 ' '' 



I d \ z ( 3 rf+2 
2F,(-r,-;-;zj+-2F, -r,-;^;z 



l + a 



1 / 3 fl' + 4 

2^1 -/", t;— ^;z 



2 + rf 



2 2 



2^- 



(17) 



Note that, since r is integer, the hypergeometric function 
2F\{—r^b\c\z) is a polynomial in z of degree r. 

In the one-dimensional case (d— 1), Eqs. (fTSl l and (fTSl l 
become 



V2r|0 =2^0 



V2r|l = 2^0 



1- 



1 + a 



a 



1-a 



2r 



2r+l 



a 



2r+l' 



(18) 



(19) 



These expressions coincide with those previously derived in 
Ref. la. 

Figure[T]displays the a-dependence of the (scaled) shifted 
collisional rates (Oj,.!,, = fi^r|.s/vo withi = 0, 1,2 and2r + 5 < 
10 for the three-dimensional case (d = 3). Of course, the 
null collisional rates Oqio = Woh = 00210 = are not plot- 
ted. Several comments are in order. Firstly, the degeneracy 
(02r~2\i = ®27|o present in the elastic limit M28II29I is bro- 
ken, yielding fi>2,_2|i < ®2fio- Analogously, the linear rela- 
tionship iia)2,| 1 ~{d— l)c02,._2|2 + ^r|o for elastic Maxwell 
particles no longer holds if a < 1, except in the case r = 1, 
where one has d(02\\ = {d—\)(i)()\2 for any a ll24l . Secondly, 
we observe that all the shifted collisional rates monoton- 
ically decrease with increasing dissipation, eventually be- 
coming negative, except those corresponding to 2r + 5 < 5. 
The physical implications of this change of sign will be dis- 
cussed in the next section. A further observation that can be 
extracted from Fig. [T]is that the impact of a on ©2^15 be- 
comes generally more pronounced as the degree 2r + 5 in- 
creases. In the case of the unshifted collisional rates V2r\s, 
a graph similar to Fig. [T](not reported here) shows a non- 
monotonic dependence on a; they first increase with in- 
creasing inelasticity, reach a maximum, and then decrease 
smoothly. In contrast to the shifted collisional rates G)2rU, 
the collisional rates V2, ly are always positive, as expected on 
physical grounds. 
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Fig. 1 Plot of (from bottom to top at a = 0) (Bj^qiq, (O^Q' ^^W ^W' 
COg,,, cOgij, COgiQ, C0^|2 tt)4i[, tt>4i(), (»2|l*' 1^212' ^"'^ ''^()I2' ^^ dimension- 
ality is rf = 3. 



4 Diverging moments in the HCS 

The Boltzmann equation for the HCS is given by Eq. ([T]i 
with V — > 0. It is more convenient to rewrite it in terms of 
the scaled distribution 

r{c{t),t) = -[2T(t)/m]'/^f{y,t), c{t) ^y/^lfW^. 
n 

(20) 
The resulting Boltzmann equation is 

5,r(c,T)+^|-.[cr(c,T)]=y*[cir,r], (21) 

2 dc 

where dr = Vod/, vZq = V2|o/vo is the reduced cooling rate, 
and J* is the dimensionless Boltzmann collision operator. 
From Eq. (ISTT i. and taking into account Eq. (|7]i, one gets the 
time evolution equation of the moments: 



<5tA/2,|.v = -«2*,-|,^2r|.v + — L ^,'r"\s's"sM*2^\s'M*2^i\^, . 

(22) 



^0 r'r",s',s" 



If the distribution function is isotropic, i.e., /*(c, t) = 
/*(c, t), then the only non-vanishing moments are M^^t^^^ix). 
We will refer to them as the isotropic moments. On the other 
hand, if the initial distribution function /* (c , 0) is not isotropic, 
the other moments, in particular M2^|.(t) and M^^^^Xz), are 



not necessarily zero. We will call anisotropic odd moments 
to M2^j.{t) and anisotropic even moments to A/jj-n/t). 

Since the time evolution of the scaled velocity moments 
in the HCS is governed by the shifted coUisional rates (02r\s, 
the fact that the latter can become negative (for a smaller 
than a certain threshold value depending on r and s) implies 
that the associated moments diverge in time. 

Among the (scaled) moments M^^,^, M^^,., and ^21;-, 
Fig. [T] shows that the lowest-degree diverging moments are 
(in the three-dimensional case) the sixth-degree moments 
M*| . . and M*|„, which diverge for a < 0.020 and a < 0. 145, 

A\ij 6|0' '=' — — ' 

respectively. Moments of higher degree diverge for smaller 
inelasticities. More specifically, M*. ., M*. . ., M*|„, M*. ., M*. . ., 

i^ -" o\i' 6\ij' 8|0' 8|(' a\ij' 

and M*o|o diverge for a smaller than 0.261, 0.331, 0.386, 
0.444, 0.482, and 0.514, respectively. In general, the larger 
the degree the larger the threshold value of the coefficient 
of restitution below which the moment diverges. Given a 
degree 2r, the isotropic moment M2^,q diverges earlier (i.e., 
with a larger threshold value a = oc2r|o) than the anisotropic 
(even) moment Ml 



2r-2\ij 



The threshold value of OI^m can be 
obtained as the solution of the equation ©2,-10 = 0. From Eq. 
( fTsT i, this is equivalent to 



^0 



a-) = l 



1 



a 






2F1 



1 d 



(23) 



Given an integer value of r, Eq. (I23I I is an equation of degree 
2r in a. 

The Boltzmann equation (I2TI) for the scaled distribution 
function /*(c,t) admits a stationary and isotropic solution 
(j)oo{c). This corresponds to a self-similar solution to the orig- 
inal Boltzmann equation where all the velocity and time de- 
pendence is encapsulated in the scaled velocity c. About ten 
years ago, Ernst and Brito 112 11122 1 conjectured that the gen- 
eral solution of Eq. ( I2TI 1 asymptotically tends to (poo{c) for 
long times. Let us loosely express this conjecture as 



limf*ic,T)^Mc), 



(24) 



where the precise meaning of the limit needs to be fixed in 
a rigorous mathematical sense. The existence of the self- 
similar solution and the convergence rate for the general 
approach to this state was first addressed in Ref. fS). How- 
ever, in that work the authors imposed conditions that were 
proven to be unnecessary in Refs. ll9l [T0l[T2l . More recently, 
proofs of the strong convergence in Sobolev and L' norms 
for small ifTTl and finite [23l| inelasticity have been pub- 
lished. Those proofs hold for any initial data (probability 
densities with bounded second-degree moments), regardless 
of being isotropic or not, but they do not imply that any 
moment of degree higher than two converges to the corre- 
sponding moment of the self-similar solution. This stronger 
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moment-to-moment interpretation of the Ernst-Brito con- 
jecture would read 



limM*,|,(T) 



dcF2r|.v(c)0oo(c) 



(25) 



As discussed below, this stronger notion of the convergence 
statement (|24] | does not hold. 

Although the explicit form of (p„{c) is not known, ex- 
cept in the one-dimensional case |l2l, it is known that it pos- 
sesses an algebraic high-velocity tail of the form (/)oo(c) ~ 
(;-d-Yom^ where 70(0;) obeys a transcendental equation ||4j 
I9lll0il2ni22ll26l . As a consequence, the isotropic moments 
^2r\o ^^^^ 2r > 70(0:) diverge. According to the strong con- 
vergence property dZSl l. this would imply that, if M2^|q(0) = 
finite, then UmT^ooM*^|Q(T) =°°if2r> 70(a). This is fully 
consistent with the fact that O2r|o < 0, so that M|^|q(t) in- 
deed diverges in time if a < a2,-\o- In fact, formally setting 
2^ = 7o in Eq. ( l23T l one recovers the transcendental equation 
for 7o derived by an independent method Il4l l2ni22ll26l . 

The interesting point is that, as shown above, the anisotropic 
moments Mi , . and Mi , . . can also diverge, unless they are 
zero in the initial state. The possibility that limT^ooMj^ , .(t) = 
°o and limT^ooMj^j. .(t) = °° contradicts the strict moment- 
to-moment limit dZSl ), since all the anisotropic moments of 
0oo(c) vanish. Let us elaborate this result in more detail. 

In principle, we have derived Eqs. (fT5Tl-(fT7Ti for r = integer. 
However, since the hypergeometric function and the Pochham- 
mer symbols are well defined for r 7^ integer, we speculate 
that an analytic continuation of Eqs. (fTSt-lfTTJi to r 7^ integer 
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Fig. 2 Plot of (from bottom to top) 70(0:), 71 (a), and 72(a). The inset 
shows (1 — 0:^)7, versus a. The dimensionality is d = 3. 



(26) 



is possible. It is then tempting to interpret (O^iq, 0^,^-111' ^^'^ 
(0^_2\2 as the quantities governing the asymptotic time evo- 
lution of the averages M*|„ = (c^), M, 



Since «^_212 > ^k-i\i ^ ^klo^ ^^ '■^™^ °^^ *^'- 



^k-l 



Ml 



— l^k-2 I 



k\Q 



k-l\i 



Ci 



and Hm 



M] 



k-lU 



^_2|,; = \'- \i,ii-]~d ^c'^5;;)),respectively,evenif A;/2 7^ 
integer and (^ — 1 )/2 7^ integer, although a formal proof of 
this expectation is beyond the scope of this paper As said 
before, M^q ^- 00 if ^ > 70(a), where (O* .q = 0. Analo- 
gously, we can expect that the anisotropic quantities M^ j,. 
and M^_ 2 1,- diverge if A: > 71(a) andA:> 72(a), respectively, 
where 71 and 72 are the solutions to the equations o* j 1 j =0 



M\ 



0, 



lim ^rr— ^ = 0. 



A'lO 



M\ 



(27) 



k\a 



and (0, 



yi-2\2 ■ 



0. 



The functions 7o(a), 7i(a), and 72(a) are displayed in 
Fig. 12] for d ~1>. In the elastic limit a ^^ \, the three ex- 
ponents diverge as js ~ Ad /{I — a^) B1 l26l . as shown in the 
inset of Fig. 121 We observe that 70 (a) < 71(a) < 72(a). This 
implies that, at a given value of a the isotropic average M^q 
starts to diverge before the anisotropic (odd) average M^ ,,. 
does, and the latter does it before the anisotropic (even) aver- 
age M^_2|,- does. Stated differently, if we focus on the ratios 
between the anisotropic and the isotropic averages, we can 
expect the asymptotic behaviors 



M\ 



/t-lh 



M\ 



k-2\ij 



-H-2 



k-l\2 



M: 



k\0 



M: 



k\0 



Therefore, the anisotropic moments, relative to the isotropic 
moments of the same degree, asymptotically go to zero (the 
anisotropic even moments more rapidly than the anisotropic 
odd ones). From that point of view, Eq. dZ/l ) can be seen as 
a weak validation of a moment-to-moment interpretation of 
Eq. (I24I 1 for initial anisotropic distributions. 

The one-dimensional system deserves some separate com- 
ments. In that case, the self-similar solution is ^oo(c) = [2^'^ /k){1 - 
2c^y^ 121, so that 7o = 3 and the moments (c'') with k > 
3 diverge. This agrees with Eq. ( fTsT i, according to which 
ft^m — ^ ^'^^ k>3. Analogously, from Eq. ( fT9l l one gets 71 = 
3. In particular, the isotropic moment (c^) diverges, while 
the anisotropic moment (c^Cx) (proportional to the heat flux) 
keeps its initial value If3l l24l . Therefore, {c^Cx)/{c^) -^ 0. On 
the other hand, since (oZq = fi^^_i|i < for ^ > 3, there ex- 
ist two possible scenarios for the ratios (c*^^'cv)/(c*^): either 
they tend to constant values or they decay more slowly than 
exponentially. A deeper investigation is needed to elucidate 
between these two possibilities. 
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5 Conclusion 

To summarize, we have shown that the strong notion dZST l 
of the Ernst-Brito conjecture cannot be strictly true since it 
does not hold for anisotropic initial conditions. However, we 
conjecture that -^2r|s^('^)/''^2r+v|o('^) ^ ^ when s ^Q, even 
if A^2rls(^) ^" °°' ^^ shown by Eq. dZTJ i for s = i {s ~ 1) and 
s = ij {s — 2). In order to elaborate further this conjecture, 
let us decompose /* (c, t) into its isotropic, anisotropic sym- 
metric, and antisymmetric parts: 



r(c,T)=</)(c,T)+/;(c,T)+/:(c,T), 



where 



(28) 



/;(c,t)=/;(c,t)-</)(c,t), (/)(c,t) = — ydc/+(c,T), 

(29) 



/i(c,T) = l[r(c,T)±r(-c,T)]. (30) 

As a consequence, the velocity moments M^|q(t), M^_j|.(t), 

and M^_2|,- ■ ('^) ^^ related to (c, t), /* (c, t), and f\_ (c, t), 
respectively. If the "sizes" of these three contributions are 
measured through those three classes of moments, we can 
say that, as time progresses, the two anisotropic parts of /* 
become negligible versus the isotropic part, i.e., \fl (c, t) | ^ 
0(c, t) and |/^(c,t)| < 0(c,t), in the sense of Eq. (|27] |. 
Moreover, limT^c»0(c) = 0oo(c). We further speculate that 
the high-velocity tails of the anisotropic contributions tend 
to the forms 

/:(c,T) -> x_(c)c-'-«(«), /;(c,T) ^ z+(c)c-''-'^-(«\ 

(31) 

where the angular functions x^ (c) and %+ (c) depend on the 
initial conditions. A confirmation of the above expectations 
requires a more refined analysis. 
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